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1.  INTRODUCTION 

Since  a real  or  potential  cost  is  involved,  any 
item  sold  with  a warranty  must  necessarily  be 
priced  higher  than  if  it  were  sold  without  a war- 
ranty. How  much  more  the  seller  should  charge 
and  how  much  more  the  buyer  should  be  willing  to 
pay  depends  upon  the  structure  of  the  warranty 
and  the  life  distribution  of  the  item.  An  anal- 
ysis of  pro  rata  and  free- replacement  warranties 
from  both  buyer's  and  seller's  points  of  view  is 
given  by  Blischke  and  Scheuer  (1975  and  1977). 

In  this  paper  we  shall  consider  only  the  free- 
replacement  warranty  and  shall  be  concerned  main- 
ly with  the  seller's  (or  supplier's,  manufac- 
turer's, and  so  forth)  point  of  view.  Of  primary 
importance  from  this  point  of  view  is  the  long- 
run  profitability  of  the  item. 

An  important  consideration  in  analyzing  long-run 
profits  for  items  sold  under  a free-replacement 
warranty  is  the  expected  income  over  the  life 
cycle  of  the  item.  This,  of  course,  is  a function 
of  the  expected  number  of  replacement  items  sold 
over  the  life  cycle.  (We  assume  here  that  the 
buyer  purchases  an  identical  replacement  when  the 
item  in  service  at  the  end  of  the  warranty  period 
fails  and,  for  simplicity,  that  the  purchase  and 
initiation  of  operation  of  the  replacement  are 
instantaneous.)  This  expected  number,  called  the 
"renewal  function"  for  the  associated  random 
variable,  is  the  subject  of  this  investigation. 

In  the  ensuing,  we  shall  discuss  in  more  detail 
the  nature  of  the  free-replacement  warranty  and 
its  associated  costs/profits,  the  role  of  renewal 
theory  in  analyzing  warranty  policies,  and  the 
specific  renewal  function  encountered  in  the  con- 
text just  described. 

The  form  of  a renewal  function  depends  ultimately 
on  the  underlying  life  distribution  of  the  items 
in  question.  Typically  in  dealing  with  renewal 
functions,  closed  form  expressions  are  available 
only  for  a few  special  cases,  although  limiting 
results  are  quite  generally  available.  We  shall 
find  this  to  be  true  of  the  special  renewal 
function  under  consideration  here  as  well.  Ana- 
lytical results  will  be  given  for  the  exponential 
distribution  and  (although  admittedly  of  limited 
interest  as  a distribution  of  lifetimes)  for  the 
uniform  distribution.  Some  preliminary  results  of 
a numerical  investigation  of  the  special  renewal 
function  for  gamma  and  Weibull  distributed  life- 
times will  also  be  discussed. 

2.  THE  SPECIAL  RENEWAL  FUNCTION  AND  ITS  ROLE  IN 
THE  ANALYSIS  OF  WARRANTY  POLICIES 

2,1  The  Analysis  of  Warranty  Policies.  In  the 
analysis  o^  warranty  policies  given  by  Blischke 
and  Scheuer  (1975  and  1977)  the  basic  considera- 
tions were  the  comparison  of  cost  to  the  consumer. 


and  of  profit  to  the  supplier,  of  warrantied 
versus  unwarrantied  items.  In  the  present  paper 
we  shall  limit  consideration  to  the  point  of  view 
of  the  supplier.  From  his  point  of  view,  the 
cost  comparison  leads  to  the  establishment  of  a 
differential  pricing  structure  which  will  equate 
expected  long-run  profit  in  the  two  situations. 

Profit,  of  course,  is  a function  of  cost  and  in- 
come. In  our  previous  work  (Blischke  and  Scheuer 
(1975))  we  derived  the  expected  profit  per  war- 
ranty cycle.  Here  we  are  concerned  witn  the 
long-run  profit  over  the  life  cycle  of  the  item. 

This  can  be  approximated  for  relatively  long  life 
cycles  by  pursuing  an  analysis  along  the  lines  of 
our  1975  paper.  (See  especially  Sections  2.1.1 
and  2.2).  Our  present  objective  is  to  obtain  an 
exact  expression  fan  this  quantity.  A result  of 
this  type  will  also  provide  a basis  for  evalu- 
ation of  the  adequacy  of  the  approximation. 

2.1.1  The  Free-Replacement  Warranty.  The  speci- 
fic warranty  policy  under  consideration  here  is 
the  free-replacement  policy.  Under  a warranty  of 
this  type  the  supplier  provides  replacements  for 
failed  items  free  of  charge  until  a specified 
period  of  service,  W,  is  attained.  His  income 
during  this  period  is  the  price,  C,  charged  for 
the  initial  item.  His  expected  cost  is  the  sum 
of  the  cost  of  supplying  the  initial  item  and  the 
expected  cost  of  all  replacements  required  to 
provide  the  total  warrantied  service  time,  W.  In 
the  sequel  we  shall  express  this  expected  cost, 
following  Blischke  and  Scheuer  (1975),  as 

g(l  + Mj,(W)l,  where  g is  the  cost  per  unit,  X is 
the  random  1 ifetime  of  an  individual  item  and 
M„(w)  is  the  associated  renewal  function  evalu- 
ated at  W.  (In  this  expression  the  quantity 
1 + M„(W)  is  the  expected  total  number  of  items 
supplied,  that  is,  the  initial  item  plus  the  ex- 
pected number  of  replacements.) 

2.1.2  The  Excess  Random  Variable.  For  the  long- 
run  analysis  of  the  free-replacement  warranty 
policy,  it  is  important  to  note  that  no  cost  is 
incurred  and  no  income  obtained  after  W until  the 
item  in  service  at  time  W fails.  The  symbol  Y is 
used  to  denote  the  random  time  at  which  this 
event  takes  place.  This  can  also  be  expressed  as 
Y = W + y(W),  where  y(W),  the  "excess  random 
variable,"  is  the  (random)  residual  lifetime  of 
the  item  in  service  at  time  W.  This  random  vari- 
able is  key  to  the  analysis  which  follows.  It  is 
also  called  the  "excess  life"  or  "residual  life" 

(Ross  (1970),  p.44),  "remaining  life"  (Barlow 
and  Proschan  (1975),  p.  168)  and  "forward  re- 
currence time"  or  "residual  life-time"  (Cox 
(1962),  p.  27)  and  has  some  unusual  properties. 

(See,  for  example.  Feller  (1966),  Sections  1.4 
and  VI. 7.) 

2.2  The  Role  of  Renewal  Functions.  In  the  fore- 
go inq”wiTave~seenTTiat~tlie~ri?5iwal  function, 

M^(-),  of  the  basic  lifetime  random  variable,  X, 

piays  an  important  role  in  determining  expected  ;n  nrZ 
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profit  on  a per-cycle  j)asis.  In  particular,  ex- 
pected profit  per  cycle  Is  P“C-g[l+Mjj(W)] . 

We  turn  now  to  the  analysis  of  long-run  expected 
profit.  In  this  case  we  look  at  repetitions  of 
the  warranty  cycles.  The  first  such  cycle  ex- 
tends from  0 to  Y.  » W + say,  the  second 

from  Y,  to  Y,,  ana  so  forth.  SchematIcaUy,  we 
have  ^ ^ 

Expected 

Cost  g[l+M^{W)]  0 g[l+M^(W)]  0 

0 W Yj  Yj+W L 

Yj(M)  Y2(W) 


work  to  approximate  Myd).  By  this  theorem, 
MyCtl/t  -*  1/E(Y)  as  ■».  (It  can  also  be  shown 
(Cf.  Smith  (1958),  equation  (1.1))  that  E(Y)  = 
u[l  + Mj,(W)j,  which  was  the  result  actually  used 
In  our  approximation.) 

3.2  Distribution  of  Y. 

3.2.1  Distribution  of  the  Excess  Random  Variable. 
Since  Y~=Tr7'YTSTr'tHi~3TstrT5UtTorro?~Os 
simply  a translation  of  the  distribution  of  the 
excess  random  variable.  Thus  the  fundamental  re- 
sult required  Is  the  distribution  of  y(W).  There 
are  several  ways  of  expressing  this  result.  All, 
of  course,  relate  back  to  the  basic  distribution 
of  X since  we  can  also  write  y(W)  as  y(W)  * 

^Nj((W)+l  ■ “• 


Income 


C ••• 


The  total  expected  profit  Is  thus  seen  to  be  P 
times  the  number  of  expected  repetitions  of  this 
process  over  the  life  cycle,  L.  This  quantity  Is 
precisely  the  renewal  function  of  the  random 
variable  Y,  evaluated  at  L.  We  call  this  the 
special  renewal  function  and  denote  It  My(-).  We 
note.  Incidentally,  that  In  our  prevlous^work  we 
approximated  this  by  the  quantity  L/E(Y).  In 
some  examples  this  approximation  was  found  to  be 
Inadequate,  which  motivated  the  present  study. 


The  remainder  of  this  paper  is  devoted  to  an 
evaluation  of  My(.).  This  function  and  M„(-)  are 
all  that  are  required  to  complete  the  cost  analy- 
sis. 

3.  ANALYTICAL  INVESTIGATION  OF  MyC-). 

3.1  General  Renewal  Theoretic  Results.  We  begin 
with  the  basic  renewal  process  Involving  a single 
warranty  cycle.  X,  Y,  y(0,  W and  L are  as  de- 
fined previously.  Let  X,,  X.,...  be  the  life- 
times of  the  individual  items  within  a warranty 
cycle.  We  assume  that  X,,  X-,-..  are  nonnegative 
random  variables  which  afe  independent  and 

distributed  with  CDF  Fy{.).  We  write 
1=1,2,...),  SnsO,,and  Ve(X).  For 
F{-),  we  define  F'"^(-)  to  be  the 
lution  of  F(-)  with  itself,  with 


Identically 
S„=z..,  X^  ( 
ally  ^'^CDf, 
n-fold  convo 


F<°>(t) 


11  t > 0 
(0  t < 0, 


N = N(t)  = number  of  replacements  required  in  the 
Interval  [0,t],  and  M(t)  = E(N(t)). 


The  distribution  of  y Is  shown  by  Ross  (1970, 
Section  3.6)  to  be  ^ 

(3.2)  Fy(w)(t)=F^(W+t)-/3  [l-Fj,(W+t-u)]  dM^(u). 

An  equivalent  expression  In  terms  of  the  corres- 
ponding densities  f and  m (when  they  exist)  Is 
given  by  Cox  (1962,  Sectijjn  5.2)  as 

(3.3)  ('y(y)(t)=f)((W+t)-^  m(W-u)fj((u+t)du. 

3.2.2  Mixture  Representation.  It  Is  of  Interest 
to  note  that  in  addition  to  these  classical  repre- 
sentations, the  distribution  of  the  excess  random 
variable  can  also  be  expressed  as  a mixture  of 
distributions  (cf.  Bllschke  (1965)  and  (1968)), 
namely 

(3.4)  F^,y,(t)  = “ P(y(W)  <t|N(W)=n)P{N(W)=n}. 

Here  the  distribution  of  N (given  In  equation 
(3.1))  Is  the  mixing  distribution  and  the  condi- 
tional distributions  of  y given  N are  the  compo- 
nents of  the  mixture.  The  difficulty  here  Is  the 
evaluation  of  these  conditional  distributions. 
Since  the  event  { n(w)  = n } is  equivalent  to  the 
event  ( S < W,  S ^ W },  the  conditional  dis- 
tributions become 

P{Y(W)<t|N(W)«n}=P{S„^j<W+t|S„<W,  S,,^j>W} 

which  can  be  expressed  as  an  Integral  over  the 
appropriate  region  of  the  bivariate  distribution 
of  S , S j.  Except  for  a few  simple  cases,  these 
Integrals  ^are  tedious  to  evaluate. 

One  property  often  encountered  In  dealing  with 
mixed  distributions  is  that  they  may  be  multi- 
modal. This  is  indeed  the  case  for  the  distribu- 
tion of  the  excess  random  variable,  a fact  that 
became  quite  apparent  In  some  of  our  computer 
simulations . 


A well-known  general  renewal -theoretic  result  is 

(3^1)  P(N(t)  = n)  = F^'’^t)  - F^'’'^^^t). 

This  provides  an  Immediate  expression  for  Mj.(t) 
in  terms  of  the  convolutions  Fi"'(t)  and  fo^ 
My(t)  In  terms  of  F'"^(t).  We^turp  next  to  the 
problem  of  determining  Fy(-)  and  Fy  M-)- 

Many  asymptotic  results  regarding  renewal  func- 
tions are  also  available.  Of  primary  interest 
here  is  the  Elementary  Renewal  Theorem  (Ross 
(1970),  p.  40),  which  was  used  In  our  previous 


Another  property  of  mixtures  of  the  type  we  are 
dealing  with  here  is  that  the  moments  of  the  mixed 
distribution  can  be  expressed  as  weighted  averages 
of  the  moments  of  the  components.  We  have  not 
pursued  this  point  but  it  would  be  of  interest 
in  some  applications.  (For  example,  on  an  indi- 
vidual basis  one  would  be  interested  in  the  con- 
ditional expected  lifetime  of  the  item  actually  in 
service  at  the  end  of  the  warranty  period.) 

3.3  Examples . 
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3-3.1  The  Exponential  Distribution.  For  the  ex- 
ponen t i al  distribution. 


(3.5) 


fx(x) 


xe 

0 


-Xx 


X ^ 
X < 


explicit  expressions  for  all  of  the  above  are 
easily  obtained.  We  use  (3.3)  to  obtain  the  den- 
sity of  the  excess  random  variable.  The  "renewal 
density"  is  m(t)  = 1/E(X)  = x.  (See  Cox  (1960), 
p.  30.1  Thus,  for  tX), 

which  is, of  course,  a well-known  result.  The 
density  of  Y is  simply  a translated  exponential. 
The  n-fold  convolution  of  this  is  a translated 
gamma,  with  CDF 

0 y<nW 

(3.6)  fI"I(,)  . ^ 

i =0  • y^nW . 

In  writing  the  renewal  function,  it  will  be  con- 
venient to  express  L as  a multiple  of  W,  say 
L = £W.  We  then  obtain,  from  equations  (3.1) 
and  (3.6), 

P(Ny(iiW)=n)  = e~^^‘^~”~^^“y:xVii-n-l)^W^ 

' i=0 


^ ! 


- e 


-x(t-n)W  "y;\i(i-n)’w’ 


i=0 


i! 


for  n = 0,l,...,i-l.  Finally,  the  special  renew- 
al function  is  found  to  be 

(3.7)  MJzW)  = (t-1) 

j=l  i*0 


i! 


3.3.2  The  Uniform  Distribution.  Although  the 
uniform  distribution  is  admittedly  of  limited  in- 
terest as  a life  distribution,  it  is  a convenient 
and  nontrivial  example  to  illustrate  the  mixture 
formulation.  (In  the  exponential  case,  all  com- 
ponents of  the  mixture  are  the  initial  distribu- 
tion of  (3.5) . ) Take 


(3.8)  fx(><)  = 


i 0 < X < 9 
0 otherwise. 


It  seems  sensible  to  assume  that  e > W since 
otherwise  replacements  are  required  with  proba- 
bility one . 

We  proceed  to  express  the  distribution  of  the 
variable  as  a mixture.  The  mixing  distribution 
is  the  distribuiton  of  N.  We  find 


P(N*0)  = P(X^>W)  = (9-W)/e  , 
P(N=1)  = P(Xj<W,  Xj+X2>W) 

’/o/w-Xj  (g2)‘'*l‘*'‘2 


29W  - W^ 


2e‘ 


and  so  forth.  The  general  result  is 


(3,9)  P(N»n)  » 

(n+Dle"^^  n!e 


The  components  of  the  mixture  are  the  conditional 
distributions  of  y(w)  given  N'n.  These  are 


(3.10)  f^(t|N=n) 


” n+1 

w"-(W-ft-e)" 

,n+l 


0 < t < e-W 


e-W  < t < e. 


eW 


n W" 


n+1 


From  equations  (3.9)  and  (3.10),  we  find  the  un 
conditional  distribution  of  the  excess  random 
variable  to  be  » i u 


(3.11) 


0<t<6-W 


^ w"-(w+t-e)"  w"  , w , 
— TK+T ^iriT'®'nTr' 


n=0  n 

n+1 


1 W/e 
9 “ 


n'e 


;;Ti?T"’-n'TT' 
e-W<t<9 
0<t<9-W 


i e“/®-  i p(W+t-9)/9 
9 9 9-W<t<e, 


from  which  it  follows  that  the  distribution  of  Y 

is 

W<y<9 


fy(y)  = 


1 W/9 
9 


i eW/8-  i 9<y<9+W. 

0 6 


The  conditional  means  of  the  excess  random  vari- 
able are  easily  seen  to  be 

-^9  - ^ ) 

E(Y(W)|N=n)  = I - ^ 5^  . 

^ 2(®  - 

One  can  use  this  result  or  equation  (3.11)  to  ob- 
tain the  unconditional  mean  of  y(W).  The  result 
is  E(y(W))  = ( 9/2)e  ' - W,  from  which  it  follows 

immediately  that  E(Y)  = ( 9/2)e"'  . 

The  convolutions  of  fy  are  rather  tedious  and  we 
have  not  pursued  the  exact  analysis.  One  could, 
however,  use  the  above  result  and  the  Elementary 
Renewal  Theorem  to  obtain  an  asymptotic  expres- 
sion for  My('). 

3.3.3  The  Gaitwa  and  Weibull  Distributions.  The 
gamma  and  Weibull  distributions,  with  respective 
densities 

(3.11) 
and 

(3.12)  fx(x) 


fj,(x) 


1 ,a-lg-x/e 

x^O 

r(a)B“ 

0 

x<0 

a ^a-lg-(x/6)“ 

x>0 

b“ 

0 

x<0, 
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are  two  of  the  more  Interesting  life  distribu- 
tions. Unfortunately  both  are  often  difficult  to 
deal  with  analytically.  This  is  certainly  the 
case  here;  general,  closed-form  expressions  for 
the  basic  renewal  functions,  to  say 

nothing  of  the  special  renewal*functions, 
exist  for  neither.  (There  is  however,  a clSsed- 
form  expression  for  the  basic  renewal  function 
for  the  gamma  distribution  if  the  shape  parameter, 
a,  is  Integer- valued.  See,  for  example,  Barlow 
and  Proschan  (1965),  page  57.)  It  follows  that 
either  asymptotic  results  or  numerical  approxi- 
mations are  required.  The  basic  renewal  function 
has  been  evaluated  numerically  and  tabulated  to 
some  extent  by  Soland  (1968).  This  may  be  used 
to  approximate  My(L)  for  relatively  large  L as 
indicated  previously. 

4.  NUMERICAL  INVESTIGATION 

4.1  Structure  of  the  Numerical  Studies.  Because 
of  the  complexity  encountered  in  the  analytical 
investigation  of  the  distribution  of  the  excess 
random  variable  and  the  evaluation  of  the  special 
renewal  function,  programs  were  written  to  pro- 
vide an  opportunity  to  investigate  the  properties 
of  both  of  these  numerically.  The  basic  life 
distributions  that  can  be  used  in  the  simulations 
with  these  programs  are  the  exponential,  ganma, 
Weibull,  uniform  and  normal.  Here  we  shall  con- 
cern ourselves  only  with  the  ganma  and  Weibull 
distributions.  Some  preliminary  results  concern- 
ing the  special  renewal  function  for  these  will 
be  discussed  below.  The  purposes  of  the  special 
renewal  program  were  to  provide  a means  of  in- 
vestigating the  approximation  My(L)/L  - 1/E(Y), 
and  to  provide  a means  of  evaluating  My(L)  when 
the  approximation  is  not  adequate. 

As  noted  previously,  E(y)  * w[l  + Mj.(W)] . Mj.(-) 
is  tabulated  by  Soland  (1968)  for  tne  ganma  *and 
Weibull  distributions  for  several  choices  of  a 
and  6 selected  so  that  in  all  cases  u°l.  We  re- 
strict consideration  to  cases  of  this  type  as 
well.  The  specific  results  which  will  be  report- 
ed are  for  o»2,3,4,5  in  combination  the  appropri- 
ate values  of  B so  that  u^l  in  each  case.  All 
combinations  of  W « 0.5,  1.0  and  1.5  with  L*5, 

10  and  15  were  used..  (This  gave  warranty  periods 
less  than,  equal  to,  and  greater  than  the  mean 
life  and  life  cycles  ranging  from  3+  to  30  times 
the  warranty  period.)-  In  each  simulation  500 
repetHions  of  the  special  renewal  process  were 
performed. 

4.2  Results.  In  each  of  the  simulations  the 
average  number  of  renewals,  say  fly(L)  was  cal- 
culated (along  with  certain  additional  relevant 
summary  statistics).  The  basic  results  for  the 
gamma  distribution  are  given  in  Table  1 and  for 
the  Weibull  distribution  in  Table  2.  In  each 
case  the  values  tabulated  are  fty(L)/L.  For  com- 
parison purposes,  values  of  1/E(Y),  based  on  So- 
land's  tables  of  My(W),  are  included  as  well. 


Table  1 

Values  of  My(L)/L  for  the 
Gamma  Distribution 


Parameters 


.707  .736  ,753 
.750  .794  .805 
.792  .830  .840 
■819  .858  .868 

W = 1.0 


.491  .526  .539 
.522  .561  .573 
.533  .578  .588 
.539  .587  .602 


.350  .394  .413 
.369  .410  .428 
.377  .424  .440 
.380  .428  .443 


Table  2 

Values  of  My(L)/L  for  the  Weibull  Distribution 


Parameters | 


W * 0.5 


a 

B 

1/E(Y) 

L 5 

10 

15 

2 

1.13 

.844 

.757 

.799 

.814 

3 

1.12 

.920 

.832 

.875 

.888 

4 

1.10 

.960 

.869 

.912 

.927 

5 

1.09 

.980 

.882 

.932 

.947 

W = 

1.0 

2 

1.13 

.616 

.532 

.574 

.593 

3 

1.12 

.653 

.568 

.612 

.625 

4 

1.10 

.668 

.575 

.622 

.641 

5 

1.09 

.675 

.581 

.633 

.646 

W « 

1.5 

2 

1.13 

.469 

.372 

.421 

.436 

3 

1.12 

.480 

.388 

.431 

.446 

4 

1.10 

.482 

.396 

.431 

.449 

5 

1.09 

.485 

.397 

.432 

.453 

In  the  simulations  we  also  calculated  the  sample 
variances  of  the  number  of  renewals.  From  these 
results  one  can  estimate  the  standard  error  of 
fty(L)/L.  The  results  ranged  from  less  than  .002 
to  .009,  with  all  standard  errors  except  those  for 
combinations  of  the  smallest  values  of  W and  L 
less  than  .005.  Given  that  the  accuracy  of  the 
computer  simulations  themselves  is  adequate,  one 
can  therefore  conclude  that  we  have  the  second 
digit  determined  to  within  one  unit  or  so  except 
for  a few  cases. 

4.3  Discussion.  It  is  important  to  note  that  the 
approximation  based  on  E(Y)  is  quite  Inaccurate: 
1/E(Y)  always  overestimates  fty(L),  with  the  dif- 
ference, of  course,  decreasing  as  L increases, 
(Thus  L/E(Y)  would  consistently  overestimate  My(L) 


which  would  lead  to  an  overestimate  of  the  expect- 
ed income  over  the  life  cycle  of  the  item.)  We 
note,  however,  that  there  is  some  consistency  in 
the  amount  of  overestimation.  The  average  amounts 
by  which  1/E(Y)  exceeds  My(L)/L  for  the  Weibull 
distribution  with  L = 5,  lo,  and  15  are  .031, 

.047,  and  .090,  respectively.  For  the  gaimta  dis- 
tribution the  corresponding  values  are  .031,  .044, 
and  .086.  In  each  case  these  values  are  slightly 
less  than  .5/L.  (Overal 1,  about  .45/L).  This 
suggests  that  1/e(y)  estimates  approximately 
[.5  + My{L)]/L.  In  other  runs  this  pattern  pre- 
vailed;’in  fact,  the  difference  appears  to  ap- 
proach .5/L  more  closely  as  L increases.  For  ex- 
ample, for  the  gamma  distribution  with  a=3  and 
6=1/3,  1/E(Y)  is  .4616.  In  a run  of  1,000  repeti- 
tions with  W=1.5  and  |,=99  we  obtained  My(99)  = 
45.213,  giving  [.5  + My(99)]/99  = .4617! 

There  remain  a number  of  unanswered  questions: 

What  is  the  appropriate  correction  for  smaller  L? 
How  is  it  related  to  a,  6,  and  W?  How  about  non- 
integer values  of  a?  Values  of  u other  than  1? 
and  so  forth.  Our  tentative  conclusion,  however, 
is  that  a better  approximation  to  Myd)  is 
L/E(Y)  - .5.  This  should  be  quite  accurate  for 
large  L and  will  be  a conservative  estimate 
(which  is  usually  desirable  when  predicting  in- 
come) for  smaller  values  of  L. 

We  note,  finally,  that  there  is  some  logic  to 
this  result.  The  quantity  1/E(Y)  should  be  the 
long-run  average  number  of  warranty  cycles  per 
unit  of  time.  This  appears  to  approximate 
[.5  + My(L)]/L.  For  these  to  agree,  the  numer- 
ator of'the  last  expression  should  be  the  ex- 
pected number  of  warranty  cycles  in  L.  This  can 
be  interpreted  in  that  way  if  we  note  that  the 
total  expected  number  of  sales  is  My(L)  + 1 
(where,  again,  the  "1"  is  the  initial  sale  and  M 
is  the  number  of  renewals),  but  adjust  this  down- 
ward by  .5  because  in  the  overall  aggregate  we 
expect  to  be  right  at  the  middle  of  a warranty 
cycle  at  time  L. 
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